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1.  Introduction 


Mobile  ad  hoc  networks  are  highly  desirable  for  military  communication  systems.  The 
nodes  in  these  networks  frequently  contain  multiple  antennas  (giving  rise  to  MIMO 
systems)  and  use  frequency-hopping  (FH)  waveforms  for  low  probability  of  intercept.  In 
mobile  ad  hoc  networks,  time-delay  estimation  (TDE)  is  often  used  to  estimate  the 
locations  of  nodes  in  the  network.  Direct  sequence  or  ultra-wideband  waveforms  are 
well-suited  for  TDE  due  to  their  large  bandwidth.  However,  FH  waveforms  are  desirable 
because  they  employ  simpler  receivers  with  large  dynamic  range  and  narrow  bandwidth 
around  each  hopping  frequency.  The  use  of  FH  waveforms  for  TDE  presents  new  challenges 
that  require  the  analysis  of  the  effects  of  fading  and  the  diversity  gained  by  integrating 
multiple  hops,  and  we  describe  new  results  in  this  report. 

The  specihc  areas  addressed  in  this  report  are  as  follows. 

a.  Study  source  signal  separation  and  channel  estimation  in  MIMO  communication 
systems  when  the  signals  are  CM  and  training  symbols  are  present.  Investigate 
constrained  Cramer-Rao  (CRB)  bounds,  develop  algorithms  that  jointly  exploit  the 
CM  signal  property  and  the  training  samples,  and  evaluate  performance  with 
computer  simulations  and  experimentally  measured  data. 

b.  Study  TDE  with  FH  communications  waveforms.  Evaluate  tradeoffs  in  hop  rate  for 
fading  channels,  and  analyze  the  value  of  training  symbols  and  the  impact  of  carrier 
uncertainty. 

c.  Investigate  the  design  of  waveforms  for  quasi-synchronous  superposition  that 
maximize  coherent  gain,  for  potential  use  in  cooperative  transmission  in  ad  hoc 
networks. 

The  report  is  organized  as  follows.  Section  2  of  this  report  contains  a  review  of  performance 
bounds  for  TDE,  including  a  discussion  of  current  trends.  Sections  3  through  6  contain  an 
analysis  of  TDE  performance  bounds  and  estimators  for  signals  that  are  received  on 
parallel,  flat,  Rayleigh  fading  channels.  The  model  is  applicable  to  FH  waveforms,  where 
the  signals  are  orthogonal  in  time,  as  well  as  a  sum-of-tones  waveform  in  which  the  signals 
are  orthogonal  in  frequency.  The  maximum-likelihood  estimator  (MLE),  CRB,  and 
Ziv-Zakai  bound  (ZZB)  are  presented  for  Eve  models  of  the  channel  gain  parameters: 

(1)  known  values,  (2)  unknown,  deterministic  parameters,  (3)  zero-mean,  complex, 
Gaussian  random  variables  (Rayleigh  fading),  and  (lA),  (2A)  which  are  cases  (1)  and  (2) 
averaged  over  the  Rayleigh  fading  model  in  (3).  Results  from  Monte  Carlo  simulations  are 
presented  that  compare  the  mean-squared  error  (MSE)  of  the  MLEs  to  the  appropriate 
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bounds.  The  simulations  indicate  that  the  ZZB  provides  a  signihcantly  tighter  bound  than 
the  CRB,  and  the  ZZB  exhibits  a  signal-to-noise  ratio  (SNR)-threshold  region  that  agrees 
with  the  MSE  performance  of  the  MLE.  The  examples  are  restricted  to  a  complex  envelope 
signal  model  at  baseband  and  Rayleigh  fading.  However,  the  analysis  is  general  and  may 
be  applied  to  passband  signals  as  well  as  Rician  fading.  In  addition,  the  ZZB  expressions 
may  enable  analytical  characterization  of  the  threshold  performance  of  the  MLE. 


2.  A  Survey  of  Time  Delay  Estimation  Performance  Bounds 


In  this  section  we  review  performance  bounds,  as  well  as  some  current  trends,  in  TDE. 
Research  over  several  decades  reveals  that  a  few  key  parameters  determine  TDE 
performance.  The  most  basic  are  the  SNR,  and  the  signal  time-bandwidth  (TB)  product; 
larger  values  for  each  are  desirable.  The  CRB  reveals  asymptotic  MLE  behavior  with 
respect  to  TB  and  SNR.  At  moderate  to  lower  SNR,  TDEs  generally  break  down  as 
ambiguities  arise  due  to  increased  noise  and  the  cross-correlation  of  the  signal,  causing  the 
TDE  to  deviate  (often  quite  sharply)  away  from  the  CRB.  Because  it  is  a  local  bound,  the 
CRB  does  not  indicate  the  threshold  behavior,  and  Ziv-Zakai  and  other  bounds  have  been 
developed  to  handle  this.  When  TD  is  measured  between  multiple  sensors,  the  coherence 
between  them  can  fundamentally  limit  the  result,  an  effect  that  occurs  in  acoustics  due  to 
the  turbulent  atmosphere.  We  discuss  modihcations  to  the  classical  bounds  that 
accommodate  the  coherence  loss,  and  reveal  a  threshold  coherence  phenomenon.  When 
communications  and  other  signals  are  used  for  TDE,  they  may  have  signihcant  nuisance 
parameters,  including  carrier  uncertainty,  unknown  symbols,  as  well  as  effects  due  to  an 
unknown  channel.  New  TDE  performance  limits  reveal  the  effect  of  these  parameters  for 
various  signal  models,  including  the  impact  of  diversity  channels  on  TDE,  and  this  is 
considered  beginning  in  section  3. 

TDE  is  basic  to  many  scenarios,  being  intimately  linked  to  detection,  array  processing, 
surveillance,  synchronization  in  communications,  range  Ending  in  RADAR,  as  well  as 
geolocation  and  tracking  in  sensor  networks.  Because  of  this,  TDE  has  a  long  history.  In 
this  report,  we  survey  TDE  performance  analysis  results,  with  the  caveat  that  such  a  large 
topic  could  fill  many  more  pages,  and  many  more  papers  could  easily  have  been  cited. 
While  the  topic  is  old,  there  continues  to  be  progress  and  unsolved  problems  become 
important,  particularly  in  the  areas  of  sensor  networking,  networked  geolocation,  and 
communications.  In  addition  to  some  now  classical  results,  we  also  discuss  some  more 
recent  problems  and  partial  performance  analysis  solutions.  In  the  rest  of  this  introduction, 
we  discuss  the  relevant  parameters  and  issues  associated  with  TDE  performance. 
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When  the  signals  are  coherent  (e.g.,  identical  up  to  a  delay  on  two  or  more  sensors,  or 
employing  a  single  sensor  when  the  signal  is  known  to  the  receiver),  CRB  analysis  on  the 
delay  parameter  indicates  that  the  performance  is  fundamentally  limited  by  the  signal  TB 
product  and  the  SNR,  e.g.,  see  Whalen  {53,  31).  The  CRB  on  TDE  is  the  subject  of  section 
2.1.  Thus,  for  stationary  signals  in  a  stationary  channel,  TDE  performance  is  enhanced  by 
longer  duration,  wider  bandwidth,  and  higher  energy  or  power  signals.  Of  course,  in  most 
scenarios  the  TB,  SNR,  and  coherence  times  are  prescribed  within  some  observable  range, 
and  the  MLE  of  TD  does  not  always  achieve  the  performance  predicted  by  the  CRB. 

The  CRB  provides  a  local  bound  on  unbiased  estimators  based  on  the  curvature  of  the 
likelihood  around  the  true  parameter  value,  but  it  does  not  account  for  large  errors  that 
occur  with  ambiguity-prone  signals.  For  example,  the  TDE  is  often  obtained  from  the  peak 
of  a  correlation  operation  (this  is  the  MLE  in  many  scenarios).  The  correlation  function  is 
oscillatory  for  many  signals  of  interest,  so  a  neighboring  peak  may  be  selected  incorrectly, 
leading  to  a  large  error  in  the  TDE.  The  looseness  of  the  CRB  thus  motivated  a  search  for 
tighter  bounds  and  a  study  of  the  conditions  under  which  the  MLE  achieves  the  CRB.  The 
study  of  tighter  bounds  that  predict  TDE  threshold  behavior  is  the  subject  of  section  2.2. 
We  focus  on  the  Ziv-Zakai  and  Weiss- Weinstein  bounds  (WWBs).  Generally,  these  bounds 
may  incorporate  a  prior  on  the  TD  parameter,  they  bound  the  MSE  and  are  not  restricted 
to  unbiased  estimation,  and  are  tighter  than  the  associated  CRB  and  subsume  it.  Finding 
these  bounds  is  less  of  a  turn-the-crank  procedure  than  the  CRB. 

In  some  cases,  the  signals  arriving  at  two  or  more  sensors  may  be  only  partially  coherent. 
This  is  the  subject  of  section  2.3.  This  situation  arises,  for  example,  in  aeroacoustics  due  to 
random  fluctuations  in  the  propagation  medium.  Although  this  effect  has  long  been 
recognized  {30),  TDE  performance  effects  have  only  recently  been  fully  characterized,  by 
extending  Ziv-Zakai  analysis  to  incorporate  spatial  coherence  {25).  Here,  a  threshold 
coherence  effect  is  observed,  such  that  for  a  given  TB,  the  coherence  must  be  above  a 
threshold  to  obtain  useful  TDEs,  independent  of  SNR.  The  model  used  here  has 
connections  to  multiplicative  noise  models  used  in  electromagnetic  signal  processing. 

In  communications,  the  impact  of  TDE  is  reflected  in  the  acquisition  and  synchronization 
steps  early  in  the  receiver  processing  chain.  It  is  also  of  interest  to  exploit  communications 
devices  and  waveforms  for  TDE,  for  example,  to  perform  geolocation  in  a  network.  This 
leads  to  the  topic  of  TDE  when  diversity  channels  are  available,  and  this  is  treated 
beginning  in  section  3. 
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2.1  Cramer- Rao  Bound 


Deterministic  signal  model:  Many  early  TDE  performance  stndies  concentrated  on  the 
CRB.  Consider  the  signal  model 


r(t)  =  sit]  9)  +  n(t),  0<t<T.  (1) 

Here,  the  real-valned  deterministic  signal  s{t)  is  known  np  to  some  parameters  in  6.  The 
noise  is  assnmed  to  be  additive  white  Gaussian  (AWGN),  with  two-sided  power  spectral 
density  (PSD)  No.  The  simplest  case  is  s(t;  6)  =  s{t  —  t),  9  =  r,  and  it  is  assumed  that  the 
delay  r  is  such  that  s{t-,9)  is  completely  contained  in  the  interval  [0,T],  with  s{t)  non-zero 
in  [0,Ts].  This  continuous  time  model  ignores  any  sampling  or  quantization  error  that 
arises  in  a  discrete-time  version  of  equation  1. 


The  CRB  is  given  by  (e.g.,  see  Helstrom  and  references  therein  {15)) 


CRB(r) 


1 

SNR-/32- 


with  SNR  and  signal  energy  given  by 


SNR 


E 


E  =  [  s‘^{t)dt. 

Jo 

The  factor  (5'^  is  the  mean  square  bandwidth  of  s(t),  given  by 


(2) 

(3) 

(4) 


^  j^^s\t)dt  !Z\S{f)\M  ’ 

where  S{f)  is  the  Fourier  transform  of  s{t).  For  connection  with  the  discrete-time  case  of 
equation  1,  see  {21).  The  CRB  equation  2  is  lowered  by  increasing  the  SNR,  or  increasing 
the  signal  bandwidth.  The  result  is  not  strongly  dependent  on  the  particular  shape  of  s(t), 
but  rather  its  overall  energy  and  bandwidth. 


For  equation  1,  the  MLF  results  from  maximizing  the  likelihood  function.  The  TD  MLF 
for  this  case  is  to  choose  the  peak  output  of  the  matched  hlter,  e.g.,  {15).  More  generally, 
when  there  are  nuisance  parameters,  so  that  6^  is  a  vector  including  r  and  other  unknowns, 
then  the  matched  hlter  generalizes  to  a  search  over  possible  values  for  9  in  the  correlation 
of  the  received  signal  with  s{t]  9).  For  example,  when  the  signal  has  unknown  frequency 
shift,  perhaps  due  to  Doppler  effects,  then  the  resulting  correlation  is  referred  to  as  the 
ambiguity  function. 
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The  narrow  band  case  of  equation  1  allows  us  to  write  the  signal  as 

s(t;  6)  =  A-Re  {f{t  -  r)  exp{juj{t  -  r))}  ,  (6) 

where  f{t)  is  the  complex  envelope,  u  is  the  carrier  frequency,  and  A  is  the  amplitude.  The 
problem  of  detection  and  time-delay  estimation  with  equation  6  has  been  studied 
extensively  in  the  context  of  radar  and  communications,  e.g.,  see  {15). 

Gaussian  signal  model:  When  the  signal  model  is  not  known  precisely,  an  alternative  is 
to  assume  it  is  random  and  Gaussian.  This  model  has  been  studied  extensively,  e.g.,  in  the 
context  of  acoustics,  and  many  works  in  {6)  address  this  case.  With  Gaussian  signal  and 
noise,  second-order  statistics  are  sufficient. 


Gonsider  the  two-sensor  case,  given  by 

ri{t)  =  s{t)  +  ni{t) 

r2{t)  =  s{t  -  t)  +  n2{t).  (7) 

Signal  s{t)  is  a  realization  of  a  random  Gaussian  process,  observed  on  an  array,  with 
independent  AWGN  on  each  element.  Now,  TDE  corresponds  to  the  delay  observed 
between  the  sensors,  so  is  sometimes  referred  to  as  time  difference  of  arrival  (TDOA) 
estimation.  The  signal  is  observed  over  a  time  interval  of  length  T.  This  model  allows  us  to 
fix  the  signal  power  spectrum  while  varying  the  observation  time. 


The  GRB  on  TDE  for  this  model  is  developed  in  Ghow  and  Schultheiss  {9),  which  we 
review  for  the  case  of  two  sensors.*  Let  the  PSDs  of  the  signal  and  noise  processes  in 
equation  7  be  denoted  by  >S'(a;),  A^i(ci;),  N2{uj).  The  Fisher  information  for  TDE  with  large 
time-bandwidth  product  is,  in  general  (h), 

rji  pCSD 

JcR  {t)  =  TT  /  SNR(a;)  du,  (8) 


where  SNR(a;)  is  dehned  as 

SNR(o;)  = 


(S(u)/N,(u))(S(u)/N2(u)) 

1  +  (^(o;)/iVi(o.))  +  (,S(o;)/iV2(o^)) 


and  the  GRB  =  Gonsider  the  case  where  the  signal  and  noise  PSDs  are  flat, 

n/  N  _  f  “S',  |a;|  G  [cUc  —  5/2,  lOc  +  5/2] 

^  1  0,  otherwise  ’ 


and  A^i(ci;)  =  N2{uj)  =  No,  where  the  center  frequency  is  Uc  and  the  bandwidth  is 
5  rad/sec.  Then  the  SNR  in  equation  9  is 


SNR 


{S/No)^ 

1  +  2  (S/No) 


*A  general  expression  with  a  linear  array  is  given  in  (P). 


(9) 


(10) 


(11) 
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and  the  Fisher  information  in  equation  8  evaluates  to 


JcR  (t) 


SNR. 


(12) 


The  analysis  employs  Fourier  series,  and  T  ^  2tt/B  is  needed  to  satisfy  an  assumption 
that  the  Fourier  series  coefficients  are  independent. 1  As  the  observation  time  doubles,  the 
CRB  shifts  downward  by  3  dB.  Also,  the  CRB  scales  with  the  inverse  of  SNR  S/No  and 
with  the  inverse  of  center  frequency  uJ/. 


For  this  model  equation  7,  the  MLE  of  r  corresponds  to  choosing  the  peak  output  of  the 
cross-correlation  between  ri(t)  and  r^it),  e.g.,  see  6),  and  references  therein.  This 
theory  has  been  developed  extensively  (h),  including  estimators  that  account  for 
intervening  linear  channels,  leading  to  the  so-called  generalized  cross-correlator  (GCC). 
The  GCC  generally  assumes  known  statistics,  and  is  a  weighted  (non-parametric) 
cross-correlation  TDE,  e.g.,  see  Hero  and  references  therein  {16).  The  case  of  unknown 
statistics  was  considered  by  Friedlander  and  Porat  {11). 


2.2  Ziv-Zakai  and  Weiss- Weinstein  Bounds 


The  CRB  predicts  the  MLE  asymptotic  performance,  as  a  function  of  TB  and/or  SNR. 
However,  it  fails  to  predict  the  TDE  behavior  at  more  moderate  values  of  these 
parameters.  This  threshold  behavior  arises  from  ambiguities  that  are  particular  to  the 
signal  and  its  correlation.  For  example,  when  TDEs  arise  from  cross-correlation  with 
narrow  band  signals,  a  neighboring  peak  may  be  incorrectly  selected,  leading  to  a  large 
TDE  error,  e.g.,  see  lanniello  {18). 

Several  approaches  have  been  used  to  develop  bounds  tighter  than  the  CRB.  Chow  and 
Schultheiss  {9)  considered  this  problem,  employing  the  Barankin  bound  (BB),  specihcally, 
the  Chapman-Robbins  version  of  the  BB.  lanniello  investigated  the  problem  by  analyzing 
the  probability  of  large  errors  {18,  39,  20).  Ziv  and  Zakai  developed  an  approach  based  on 
hypothesis  testing  {58),  and  various  cases  and  improvements  were  subsequently  derived 

{41,  8,  3,  4,  46). 

The  ZZBs  have  been  found  to  be  tighter  than  the  BBs,  and  the  ZZBs  generally  converge  to 
the  corresponding  CRB,  so  that  they  subsume  CRB  analysis.  The  ZZB  and  its  variants 
provide  tight  bounds  for  threshold  behavior  of  TDEs,  predicting  when  the  MLE  breaks 
down  and  departs  from  the  CRB  with  good  accuracy  for  many  TDE  scenarios.  These 
bounds  generally  account  for  a  priori  limits  on  the  time  delay.  The  ZZBs  are  Bayesian 
bounds  that  treat  the  time  delay  as  a  random  parameter,  typically  uniformly  distributed  on 

'fPor  this  model,  Gaussianity  of  the  Fourier  series  coefficients  follows  immediately.  For  more  general  signal 
models,  Gaussianity  follows  via  a  central  limit  theorem  argument  with  large  T. 
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an  interval.  In  addition,  the  ZZBs  bound  MSE,  and  are  not  limited  to  unbiased  estimates 
as  is  the  CRB,  which  treats  the  time  delay  as  a  deterministic  parameter.  Given  prior  limits 
on  r,  it  can  be  expected  that  estimators  will  be  biased  when  r  approaches  the  boundaries. 

We  review  the  ZZB  beginning  with  the  model  in  equation  7,  with  observation  time  T, 

— T/2  <t<T/2.  The  signal  and  noise  are  assumed  Gaussian,  with  spectral  densities  S{uj), 
Ni{u),  and  N2{ui).  The  ZZB  begins  with  a  binary  decision  problem  based  on  some 
arbitrary  estimate  f  of  r. 

Decide  Hq  :  t  =  a  if  |f  —  a|  <  |f  —  a  —  0|, 

Decide  Hi  :  r  =  a  +  9  if  |f  —  a|  >  |f  —  a  —  6^1.  (13) 


Without  loss  of  generality,  the  two  hypothesized  delays  are  considered  to  be  equally  likely, 
so  the  probability  of  error  is  P{f  —  a  >  0/2\t  =  a}/2  +  P{f  —  a  —  9  <  —9/2\t  =  a  +  9}/2. 
Now,  the  error  in  estimating  r  by  f  is  e  =  f  —  r. 

We  are  interested  in  the  mean  square  error  =  E{ee^)  with  r  uniformly  distributed  within 
the  interval  [—D/2,  D/2],  and  9  within  [0,D].  The  basic  ZZB  is  given  by  (5) 

1  nD  pD/2—0 

e^>—  9  d9  Pe{a,a  +  9)da,  (14) 

^  Jo  J-D/2 

where  Pe{a,  a  +  9)  is  the  minimum  attainable  probability  of  error  associated  with  the 
likelihood  ratio  test  between  the  two  hypothesized  delays.  It  immediately  follows  that 

Pe{a,  a  +  9)<  i{D(e  >  9/2\t  =  a)  +  P{e  <  -9/2\t  =  a  +  9)}.  (15) 

An  improvement  is  possible  by  noting  that  De(a,  a  +  9)da  is  a  nonincreasing 

function  of  9  {3).  Therefore,  by  using  a  nonincreasing  function  V[-]  that  hlls  the  valleys  of 
the  bracketed  function,  an  improved  ZZB  is  given  by 

1  pD  pDf2—0 

e^>—  9V  /  Pe{a,a  +  9)da  d9.  (16) 

^  Jo  J-D/2 

If  Pe{a,  a  +  9)  =  Pe{9),  so  it  is  independent  of  a,  then  equation  16  may  be  simplihed  to  (43) 

^  ^  9  Vl(D  -  9)P,(9)]  d9.  (17) 

When  the  observation  time  T  is  large  compared  to  the  correlation  time  27r/D  where  B  is 
the  bandwidth,  then  Pe{9)  is  closely  approximated  by 


Pe(0) 


^a{e)+b{9) 


(18) 
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where 


a{9) 


T 


ln[l  +  SNR(cc;)  sin^(a;6'/2)]  du, 


m 


T  SNR(o^)sin2(o;0/2) 

271  Jq  1  +  SNR(a;)  sin^(a;6'/2)  ^ 


SNR(ct;)  is  defined  in  equation  9,  and 


1 

•P(x)  —  _  /  du. 


(19) 

(20) 


(21) 


The  derivation  of  equation  18  assumes  large  TB.  As  noted  earlier,  this  condition  is 
common  in  TDE  performance  analysis,  and  is  typically  met  in  practical  scenarios;  poor 
estimates  may  result  in  practice  if  this  condition  is  not  satisfied. 


The  ZZB  generally  requires  numerical  evaluation  to  approximate  the  integrals  in  equation 
19,  20,  and  17.  However,  Weiss  and  Weinstein  derived  closed-form  expressions  for  disjoint 
bound  segments,  separated  by  thresholds,  for  narrow  band  (4S)  and  wideband  (49) 
waveforms.  For  example,  when  the  signal  and  noise  PSDs  are  fiat  with  Si(u)  =  82(00)  as  in 
equation  10  and  small  fractional  bandwidth  B/uOc  1,  it  is  shown  in  (49)  that  the  CRB  is 
attainable  only  if  the  SNR  in  equation  11  exceeds  a  threshold  value  given  by 


SNRth  = 


(22) 


where  denotes  the  inverse  of  the  function  in  equation  21.  The  SNR  threshold  in 
equation  22  is  a  decreasing  function  of  TR  and  B/uOc-  The  closed-form  ZZB  expressions 
developed  in  (48,  49)  are  valuable  for  simplified  analytical  performance  study  when  the 
fractional  bandwidth  is  small  or  moderate.  Simplified  ZZB  expressions  for  large  fractional 
bandwidth  (the  ultra-wideband  case)  have  also  been  developed  (17).  The  ZZB  in  the  above 
works  is  generally  restricted  to  estimation  of  scalar  random  variables;  it  was  extended  to 
arbitrarily  distributed,  vector  random  variables  by  Bell  et  ah  in  (2). 


Simulation  examples  comparing  cross-correlation  TDEs  and  the  corresponding  ZZBs  show 
close  agreement.  For  example,  using  realizations  of  Gaussian  signals  with  fiat-topped 
spectra,  and  given  sufficient  time-bandwidth  (TB  >  50,  say),  the  ZZB  accurately  predicts 
threshold  behavior  for  narrow  band  and  wideband  cases  (19),  as  well  as  for  ultra-wideband 
(17).  These  results  suggest  that,  for  equation  7,  the  cross-correlator  is  essentially  optimal 
even  in  the  non-asymptotic  regime. 


So  far  we  have  focused  on  the  ZZB,  which  can  be  readily  evaluated  numerically,  and 
approximations  are  available  for  some  cases.  Another  bounding  approach  is  available,  due 
to  Weiss  and  Weinstein,  that  is  referred  to  as  the  WWB  (44,  50,  51,  45).  The  WWB  is  a 
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family  of  MSE  bounds  on  moments  of  random  parameter  estimators;  the  most  general 
statement  is  in  (45).  Unfortunately,  space  limits  our  ability  to  give  a  full  statement  of  the 
WWB.  The  WWB  does  not  require  regularity  conditions,  incorporates  a  priori  distribution 
on  the  parameter,  and  is  not  restricted  to  unbiased  estimates.  Application  of  the  WWB 
requires  exploration  of  test  points,  and  optimization  over  parameters  in  the  bound,  and 
this  must  be  carried  out  for  each  case  of  interest.  Simple  expressions  for  the  WWB  for  time 
delay  have  been  derived  for  square  pulses  (50),  and  narrow  band  pulses  with  a  trapezoidal 
envelope  (51).  Note  that  regularity  is  violated,  e.g.,  with  a  square  pulse,  so  that  the  CRB 
does  not  exist. 

For  the  Gaussian  case,  the  ZZB  is  readily  evalnated  nnmerically  via  the  above  eqnations, 
and,  as  we  have  noted,  some  specihc  cases  have  been  developed  for  the  ZZB  and  WWB. 
The  ZZB  and  WWB  provide  signihcantly  tighter  performance  bounds  than  the  CRB,  but 
also  require  more  effort,  whereas  the  CRB  is  more  of  a  turn-the-crank  analytical  procedure. 
For  both  the  ZZB  and  the  WWB,  it  may  be  possible  to  hnd  simplihed  expressions  that  are 
easily  manipnlated  and  provide  valuable  insight.  The  choice  of  test  points  is  important  for 
the  WWB  (and  BB),  and  this  can  make  evalnation  of  the  bonnd  signihcantly  more 
complex  as  more  parameters  are  added  to  the  problem.  The  vector  ZZB  also  faces  a 
complexity  increase  due  to  the  maximization  needed  (2).  An  interesting  comparison  of  the 
ZZB  and  WWB  for  matched  held  processing  with  mnltiple  parameters  shows  that,  for  this 
problem,  the  ZZB  provides  a  somewhat  better  characterization  of  threshold  behavior  than 
does  the  WWB  (57). 

The  WWB  applies  to  vector  parameters  and  fnnctions  of  the  parameters,  and  the  WWB 
family  was  recently  extended  by  Renaux  et  ah,  (36).  Barankin  bonnds  have  been  developed 
by  Zeira  and  Schnltheiss  (55,  56)  for  active  ranging  systems  with  mnltiple  echoes,  and 
Schnltheiss  et  al  stndied  TDE  in  the  non-Gaussian  case  in  (40). 

2.3  Impact  of  Coherence 

The  TDE  model  in  eqnation  7  assnmes  that  the  random  signal  s(t)  is  identical  at  the  two 
sensors.  As  mentioned  above,  the  model  has  been  extended  to  allow  the  signal  components 
at  the  two  sensors  to  differ  by  an  unknown,  deterministic,  linear  hltering  operation,  e.g., 
see  ( 1 6)  and  references  therein.  However,  in  some  scenarios,  random  phenomena  in  the 
propagation  medinm  cause  a  loss  in  the  coherence  (i.e.,  correlation)  between  the  signals 
that  arrive  at  a  pair  of  spatially  separated  sensors.  In  this  section,  we  review  the  impact  of 
signal  coherence  loss  on  TDE,  and  we  present  CRB  and  ZZB  resnlts  that  are  explicitly 
parameterized  by  the  signal  coherence  (25,  24).  TDE  performance  is  signihcantly  impacted 
by  imperfect  signal  coherence,  and  the  ZZB  analysis  reveals  that  there  exists  a  “threshold 
coherence”  valne  that  must  be  exceeded  in  order  for  TDE  performance  to  achieve  the  CRB. 
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Signals  propagating  in  many  media  undergo  spatial  coherence  losses.  In  aeroacoustics 
(sounds  propagating  through  air),  signal  coherence  degrades  with  increased  spatial 
separation  between  the  sensors  due  to  random  scattering  caused  by  atmospheric 
turbulence.  Physics-based  statistical  models  for  the  coherence  as  a  function  of  frequency, 
range,  weather  conditions,  and  sensor  separation  are  reviewed  and  applied  to  signal 
processing  performance  analyses  in  {54,  21,  38,  26,  28).  Scattering  and  coherence  losses 
occur  in  underwater  acoustics,  e.g.,  see  {34)  for  an  early  analysis  of  the  CRB  on  TDE  with 
scattering,  with  more  recent  work  in  {5,  22)  and  elsewhere.  In  medical  ultrasonics  (sounds 
in  the  MHz  range  propagating  through  living  tissues),  the  CRB  on  TDE  has  been  analyzed 
for  partially  correlated  speckle  signals  {43,  7).  Spatial  coherence  losses  occur  with 
electromagnetic  signals  due  to  local  scatterers  near  the  transmitter,  e.g.,  {35),  and  narrow 
band  angle-of-arrival  estimation  with  imperfect  spatial  coherence  has  been  studied. 

We  extend  the  model  in  equation  7  as  follows  to  study  TDE  with  partial  signal  coherence 
{25): 


ri{t)  =  si{t)  +  ni{t) 

r2{t)  =  S2{t -t) +  n2{t),  (23) 


where  si(t),S2(t)  are  real-valued,  continuous-time,  zero-mean,  jointly  wide-sense  stationary, 
Gaussian  random  processes  with  PSDs  S'i(a;),  5*2  (cn)  and  cross-spectral  density  (CSD) 
S'i2(ci;).  The  AWGN  processes  ni{t),  n2{t)  have  PSDs  iVi(ci;),  N2{co),  and  they  are 
independent  from  each  other  and  from  Si{t),  S2{t).  The  magnitude-squared  coherence 
function  for  the  signals  Si{t),  S2{t)  is 


l7s(o;)r 


S'i(cn)  S'2(i:n)  “ 


(24) 


Coherence  magnitude  =  1  corresponds  to  perfect  correlation  between  the  signal  components 
at  the  sensors,  so  si{t)  and  S2{t)  are  related  by  a  linear  hltering  operation.  Coherence 
magnitude  <  1  models  frequency-selective,  random  effects  in  the  propagation  paths  from 
the  source  to  the  sensors.  The  CSD  for  the  observations  ri{t),  r2{t)  in  equation  23  is 


Ri2{uj)  =  Si2{uj)  exp{juT), 


(25) 


and  the  coherence  for  ri{t),  r2{t)  is  dehned  analogously  to  equation  24  and  may  be 
expressed  as 


l7.(o^)l^ 


l7s(o;)r 


1  + 


J 


1  + 


\N2HJ 


(26) 
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TDE  performance  bounds  (CRB  and  ZZB)  are  characterized  by  the  following  SNR-like 
quantity, 


SNR(o;) 


l-|7.(o;)r 


ll7s(o;)r  L 
l7s(o^)r 
l-|7s(o;)r 


(27) 


(28) 


which  reduces  to  the  SNR  in  equation  9  when  |7s(ci;)|  =  1.  Note  that  signal  coherence  loss 
(|7s(ct;)|  <  1)  severely  limits  the  SNR  quantity  in  equation  27,  where  the  upper  bound  in 
equation  28  corresponds  to  large  SNR  per  sensor,  (S'l/iVi),  {S2/N2)  00. 


The  Fisher  information  for  TDE  with  TB  ^  271  and  partial  signal  coherence  is  obtained 
using  equation  27  in  equation  8.  For  the  case  in  which  the  signal  and  noise  PSDs  are  flat 
with  S'i(ci;)  =  S2{oj)  as  in  equation  10,  we  can  use  equation  28  to  bound  the  CRB  as 


CRB  (r)  > 


(i/h.P)  -  1 


2-?  (if) 


1  +  ^ 


(29) 


Note  that  the  CRB  is  limited  by  the  signal  coherence  and  is  not  improved  by  increased 
SNR  per  sensor,  (S'l/iVi),  {S2/N2). 


The  general  ZZB  framework  was  reviewed  in  equations  13-21,  and  next  we  apply  the  ZZB 
to  the  passband  case  in  which  the  signal  and  noise  PSDs  are  flat  with  S'i(a;)  =  5*2  (cu)  as  in 
equation  10,  7s (cu)  =  7*,  iVi(a;)  =  N2{uj)  =  No,  and  the  fractional  bandwidth  is  small, 

B/ujc  1.  Then  the  SNR  in  equation  27  is 


SNR  = 


(30) 


For  perfectly  coherent  signals,  the  ZZB  analysis  in  {^9)  shows  that  the  CRB  is  attainable 
only  if  the  SNR  in  equation  11  exceeds  the  threshold  value  in  equation  22.  The  ZZB  is 
extended  to  partially  coherent  signals  [25)  by  using  equation  30  in  the  inequality 
SNR  >  SNRth,  which  combined  with  |7sp  <  1  yields  the  following  threshold  values  for 
S/No  and  coherence  that  must  be  exceeded  for  CRB  attainability: 


> 


1  + 


1 


SNRth 


-  1 


(31) 


|77> 


1  + 


1 


> 


SNRth 


1  + 


1 


as 


SNRth 


00. 


(32) 
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THRESHOLD  COHERENCE,  I  Y  | 


Note  that  the  threshold  coherence  value  in  equation  32  is  a  function  oiTB  and  fractional 
bandwidth  (B/uc)  through  the  formula  for  SNRth  in  equation  22. 

Figure  la  contains  a  plot  of  the  threshold  coherence  in  equation  32  as  a  function  of  the  TB 
product  for  several  values  of  fractional  bandwidth  (B/uc)  with  (S/No)  cxd.I-  Note  from 
hgure  la  that  for  threshold  coherence  values  in  the  range  from  about  0.1  to  0.95,  each 
doubling  of  the  fractional  bandwidth  reduces  the  required  TB  hj  a  factor  of  10. 

Figure  lb  shows  results  from  a  simulation  example  comparing  the  root-mean-square  (RMS) 
error  of  a  cross-correlation  TDE  estimator  with  the  CRB  from  equation  12  for  several 
values  of  signal  coherence,  7s.  The  signals  are  Gaussian  random  processes  with  bandwidth 
B  =  27r30  rad/s  centered  at  Uc  =  27rl00  rad/s.  The  signal,  noise,  and  coherence  are  flat 
over  the  frequency  band,  with  S/No  =  100  (20  dB)  and  the  observation  time  is  T  =  3  s. 
The  threshold  coherence  for  this  case  is  0.41,  from  equations  32  and  22.  Note  in  figure  lb 
that  the  simulated  RMS  error  on  TDE  diverges  sharply  from  the  CRB  very  near  to  the 
threshold  coherence  value  of  0.41,  illustrating  the  accuracy  of  the  threshold  coherence  in 
equation  32.  Other  examples  of  the  threshold  coherence  effect  for  accurate  TDE  with 
measured  aeroacoustic  data  are  available  in  {25,  1). 


S  RMS  TDE  (WIDEBAND):  f^=  100  Hz,  Af  =  30  Hz 


TIME  X  BANDWIDTH  COHERENCE 


(a)  (b) 

Figure  1.  (a)  Threshold  coherence  value  from  equation  32  versus  time-bandwidth  product  TB  for 
several  values  of  fractional  bandwidth  B/ujc  with  S/N  oo.  (b)  Comparison  of  simulated 
RMS  error  for  TDE  with  CRBs  and  threshold  coherence  value  for  signals  with  bandwidth 
B/{2Tr)  =  30  Hz  centered  at  a'c/(27r)  =  100  Hz. 

Rt  is  shown  in  (25)  that  the  threshold  coherence  values  are  nearly  equivalent  for  {S/No)  >  10  dB  and 
{S/No)  oo  when  TB  27r. 
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The  generalized  cross-correlator  may  be  used  for  TDE  with  partial  signal  coherence, 


arg  max 


Ri2{!^)G{(jj)  exp  (jcu.^)  do;, 


(33) 


where  -Ri2(i^)  is  an  estimate  of  the  CSD  in  equation  25  and  G{(jj)  is  a  weighting  function. 
The  ML  estimator  uses  weighting  function  Gml{^^)  =  SNR(a;)/  |i?i2(ci;)|,  while  the  the 
Wiener  processor  {16)  may  be  more  robust  with  Gwp{i^)  =  |7r(i^)|^- 


2.4  Sensor  Networks  and  Geolocation 


TDE  has  also  received  recent  renewed  interest  in  the  context  of  sensor  networks,  as  well  as 
geolocation  of  mobile  ad  hoc  network  nodes.  Many  envisioned  applications  require  node 
localization  and/or  tracking  of  position,  and  time  delay  estimates  between  nodes  may  be 
translated  into  accurate  position  estimates.  This  motivates  the  use  of  communications 
waveforms  for  TDE  in  wireless  networks,  e.g.,  see  {33).  Here,  parametric  models  are 
applicable  that  incorporate  nuisance  parameters,  such  as  carrier  frequency  and  phase 
offsets,  unknown  symbols,  and  unknown  channels.  CRBs  for  TDE  using  these  models  are 
important  in  the  study  of  synchronization  for  communications  {42,  32).  In  the 
communications  scenario,  typically  medium  to  high  SNR  is  of  interest,  otherwise  the  SNR 
may  be  insufficient  to  provide  a  meaningful  bit  error  rate.  However,  for  TDE  based 
geolocation,  low  SNR  and  threshold  behavior  is  currently  of  interest  {29). 


3.  TDE  on  Parallel,  Flat  Fading  Channels 


The  remainder  of  this  report  presents  an  analysis  of  TDE  performance  bounds  and 
estimators  for  signals  that  are  received  on  parallel,  flat,  Rayleigh  fading  channels.  The 
model  is  applicable  to  EH  waveforms  as  in  {29),  where  the  signals  are  orthogonal  in  time, 
as  well  as  a  sum-of-tones  waveform  in  which  the  signals  are  orthogonal  in  frequency.  The 
MLE,  CRB,  and  ZZB  are  presented  for  Eve  models  of  the  channel  gain  parameters: 

(1)  known  values,  (2)  unknown,  deterministic  parameters,  (3)  zero-mean,  complex, 
Gaussian  random  variables  (Rayleigh  fading),  and  (lA),  (2A)  which  are  cases  (1)  and  (2) 
averaged  over  the  Rayleigh  fading  model  in  (3).  We  consider  cases  (lA)  and  (2A)  because 
they  are  sometimes  used  when  case  (3)  is  intractable.  However,  in  our  model,  all  cases  are 
tractable,  allowing  comparison  of  the  various  bounds.  We  present  results  from  Monte  Carlo 
simulations  that  compare  the  MSE  of  the  MLEs  to  the  appropriate  bounds.  The 
simulations  indicate  that  the  ZZB  provides  a  significantly  tighter  bound  than  the  CRB, 
and  the  ZZB  exhibits  a  SNR-threshold  region  that  agrees  with  the  MSE  performance  of  the 
MLE.  The  examples  in  this  report  are  restricted  to  a  complex  envelope  signal  model  at 
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baseband  and  Rayleigh  fading.  However,  the  analysis  is  general  and  may  be  applied  to 
passband  signals  as  well  as  Rician  fading.  In  addition,  the  ZZB  expressions  may  enable 
analytical  characterization  of  the  threshold  performance  of  the  MLE. 


4.  Signal  Model 


A  general  model  for  the  complex  envelope  for  TDE  over  N  parallel  channels  is 

ri{t)  = '^iSi{t  -  t)  +  ni{t),  i  =  l,...,iV.  (34) 

The  Si{t)  are  known  waveforms,  ni{t)  are  complex,  AWGN  processes  with  zero  mean  and 
two-sided  PSD  Mo-,  and  each  7*  is  a  complex  scalar  representing  the  gain  of  the  channel. 

We  consider  three  models  for  the  vector  of  complex  channel  gains,  7. 

(1)  7  is  known. 

(2)  7  is  an  unknown,  deterministic,  complex-valued  parameter  vector. 

(3)  Each  7*  is  a  complex,  circular,  Gaussian  random  variable  with  zero  mean  and  variance 

^  {\li\^}  —  ^1)  independent  for  i  =  1, . . . ,  iV.  This  corresponds  to  Rayleigh  fading 
that  is  independent  from  channel  to  channel,  and  is  denoted  by  7  ~  GN  (0,  I). 

The  average  signal  power,  is  assumed  known,  and  I  is  the  identity  matrix. 

(lA)  Gase  (1)  averaged  over  the  Rayleigh  fading  distribution,  so  7  ~  GN  (0,0-^/).  This 
corresponds  to  a  fading  channel  and  a  receiver  that  has  a  perfect  estimate  of  7,  so 
that  the  TDE  is  estimated  with  a  coherent  combination  of  the  diversity  channels. 

(2A)  Gase  (2)  averaged  over  the  Rayleigh  fading  distribution,  so  7  ~  GN  (0,(7^/).  This 
corresponds  to  (imperfect)  estimation  of  7  on  each  realization,  and  then  averaging 
over  the  Rayleigh  fading  distribution. 

Gase  3  is  different  from  case  2A  in  that  7  is  not  estimated  in  case  3  since  it  is  modeled  as  a 
random  vector.  It  is  also  possible  to  model  the  channel  gain  as  7*  =  6jexp(j0j)  where  hi  is 
an  unknown,  deterministic,  positive  parameter,  and  9i  is  random  with  uniform  distribution 
over  [0,  27r)  rad  and  independent  for  i  =  1, . . . ,  iV.  The  performance  of  this  case  is 
intermediate  between  cases  2  and  3,  and  we  do  not  consider  it  here. 

Our  objective  is  to  study  the  MLE,  GRB,  and  ZZB  on  the  time-delay  parameter,  r  in 
equation  34,  for  the  hve  channel  model  cases.  The  details  are  contained  in  the  succeeding 
sections,  but  the  results  are  summarized  in  table  1,  from  which  we  make  the  following 
observations. 
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Table  1.  Relations  between  the  MLEs,  CRBs,  and  ZZBs  for  the  hve  channel  models. 


Channel  Model 

1 

2 

3 

lA 

2A 

MLE 

CRB 

ZZB 

Coh  MF 
CRBi 
ZZBi 

NonCoh  MF 
CRB2  =  CRBi 
ZZB2 

NonCoh  MF 
CRB3 
ZZB3 

Coh  MF 

MCRB 

ZZBi^ 

NonCoh  MF 
ACRB  =  MCRB 
ZZB2^  =  ZZB3 

Summary  of  channel  model  cases: 

(1)  7  is  known. 

(2)  7  is  unknown,  deterministic. 

(3)  Rayleigh  fading:  7  ~  CN  (O,  l) . 

(lA)  Case  (1)  averaged  over  7  ~  CN  (O,  a'^  J). 
(2A)  Case  (2)  averaged  over  7  ~  CN  (O,  a'^  l). 


•  The  MLE  is  a  sum  of  coherent  matched  hlters  (MFs)  in  cases  1  and  lA,  since  7  is 
known.  The  MLE  is  a  sum  of  noncoherent  matched  hlters  in  cases  2,  3,  and  2A. 
Although  the  MLEs  are  identical  for  several  cases  in  table  1,  the  performance  of  the 
MLE  varies  from  case  to  case,  except  that  the  performance  is  the  same  in  cases  2A 
and  3. 

•  The  CRBs  for  cases  1,  2,  and  3  are  generally  distinct.  However,  for  the  known  signal 
waveforms  considered  here  that  have  even  symmetry,  the  r  and  7  parameters  are 
decoupled  in  the  Fisher  Information  Matrix  (FIM),  so  CRB2  =  CRBi.  The  CRB  is 
not  dehned  for  cases  lA  and  2A,  but  the  so-called  modihed  CRB  (MCRB)  (13,  12) 
and  asymptotic  CRB  (ACRB)  {37)  dehnitions  are  similar  to  these  cases.  We  show 
that  CRB3  — ACRB  =  MCRB  at  high  SNR,  and  that  the  MCRB  corresponds  to 
CRBi  =  CRB2  when  the  fixed  channel  gains  |7jp  are  replaced  by  the  ensemble 
averages,  E{|7jp}  =  a^. 

•  The  ZZB  is  unique  in  four  of  the  hve  channel  model  cases,  and  the  ZZB  is  tractable 
in  the  four  distinct  cases. 

•  The  hve  cases  can  be  unihed  with  a  Rician  fading  model,  and  the  ZZB  is  tractable  for 
this  case  as  well.  This  extension  would  be  a  useful  item  for  future  research. 

When  the  model  in  equation  34  is  applied  with  EH  waveforms  at  baseband,  then  the 
complex  envelope  of  the  narrowband  frequency  hopping  waveform  in  the  Rh  hop  is 

K 

Si{t)  =  ak,ih[t  -  {i  -  1)LT  -  kT] ,  i  =  (35) 

k=-K 


15 


where  the  pulse  shape  h{t)  is  real-valued  and  symmetric,  T  is  the  symbol  duration,  ak,i  is 
the  kih  symbol  in  the  ith  hop,  and  each  hop  has  L  =  2K  +  1  symbols  and  dwells  for  LT 
seconds.  The  symbols  are  assumed  to  be  known  and  drawn  from  a  constellation  such  as 
PSK  with  unit  magnitude,  |afc,i|  =  1.  Over  N  hops,  a  total  oi  Q  =  LN  symbols  are 
transmitted  for  a  total  signal  duration  of  QT  =  LNT  second.  The  pulse  h{t)  is  a 
square-root,  raised-cosine  pulse  with  unit  energy,  so  its  autocorrelation  is  a  Nyquist  pulse: 


hit)  hit  +  kT)  dt 


5 


k, 


where  6k  is  the  Kronecker  delta.  When  the  model  in  equation  34  is  applied  to  a 
sum-of-tones  waveform,  the  signals  in  equation  34  have  the  form 


(36) 


Siit)  =  hit)  exp{jnfit  +  0*),  i  =  (37) 

where  hit)  is  the  pulse  shape  and  fi,  (f)i  are  the  frequency  and  phase  of  the  ith  tone.  In  this 
report,  results  are  developed  as  much  as  possible  for  an  arbitrary  signal  shape  Siit),  and 
then  specialized  to  the  baseband  FH  model  in  equation  35.  The  approach  is  applicable  to 
the  sum-of-tones  model  in  equation  37,  and  it  will  be  interesting  to  carry  out  the  detailed 
analysis  for  this  case  in  future  work. 


We  conclude  this  section  with  a  statement  of  the  probability  density  function  (pdf)  for  the 
model  equation  34  for  two  cases:  the  conditional  pdf,  pc  {R  |  7;  t),  where  7  is  a  hxed 
vector,  and  the  unconditional  pdf,  pu  {R  \t)  =  E~y  {pc  {R  \  7;  r)},  that  is  averaged  over 
7  ~  CN  (0,  (jg  I).  The  autocorrelation  functions  of  the  signals  siit), . . . ,  sisfit)  are  dehned  as 


PiiO  =  /  Siit  -  Siit)  dt,  i  =  l, 


iV, 


(38) 


where  *  denotes  complex  conjugate  and  the  limits  on  all  integrals  are  (— cxd,  cx)).  The 
coherent  and  noncoherent  matched  hlters  applied  to  channel  i  are  dehned  as  follows. 


(0  =  y  hi  Siit  -  0*  riit)  dt 

y.NONCOHg)  ^  I  /■  r,(t)  dt  . 


(39) 

(40) 


Let  Vi  be  a  vector  representation  of  rj(t),  and  dehne  R  =  [ri, . . . ,  ri\f].  Then  the  conditional 
pdf  for  R  has  the  following  form,  where  factors  independent  of  7,  r  are  suppressed: 

N 


PciRlr^r)  oc  n  exp  [21t^°^(r)  -  |  • 

i=l  ■' 

The  unconditional  pdf  is  obtained  by  observing  that  7  CN  (0,  (Tg  I)  implies  that 
Vi  ~  CN  (0,  Sj),  where  the  covariance  matrix  and  inverse  are 

+  Afol 


(41) 


1 

K 


Sj  =  Si(r)s,:(r) 


= 


\H 


Mo  +  ’ 
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(42) 

(43) 


where  Sj(r)  is  the  vector  representation  of  Si(t  —  r)  and  ^  is  the  conjugate-transpose 
operation.  Then,  ignoring  terms  and  factors  that  are  independent  of  r,  the  unconditional 
pdf  is  proportional  to 


N 

(i?  I  r)  oc 

i=l 


■Afo  +  o-'^PiiO) 


_  yNONCOH 


We  make  the  following  assumptions: 


(44) 


1.  The  signal  energy  is  the  same  on  each  channel, 

Pi{0)  =  j  \si{t)f  dt  =  p{0),  i  =  (45) 

Also,  in  order  to  simplify  some  of  the  formulas,  we  will  assume  that  the  signal 
autocorrelation  functions  in  equation  38  are  the  same  on  each  channel, 

p*(0=p(0,  t  =  (46) 

2.  The  signals  have  the  property 

J  Si(t)  Si(t)*  dt  =  0,  i  =  (47) 

where  Si{t)  =  ^Si{t).  This  property  holds  for  the  sum-of-tones  model  in  equation  37 
as  is  approximately  true  for  the  FH  model  in  equation  35,  as  will  be  shown  in  item  4 
below. 

3.  The  mean-squared  signal  bandwidth  is  the  same  on  each  channel  and  is  denoted  by 
Bg  Hz,  with  dehnition 

f  ^  /  |s.(*)t  *  =  ^  df.  (48) 

where  Si{f)  is  the  Fourier  transform  of  Si{t). 

4.  For  the  FH  model  in  equation  35,  the  symbols  ak,i  are  known  and  hxed,  so  the 
bounds  (CRB  and  ZZB)  are  a  function  of  the  particular  symbol  sequence.  However, 
the  CRB  and  ZZB  expressions  are  considerably  simplihed  if  we  assume  that  time 
averages  over  the  symbols  are  equal  to  the  corresponding  ensemble  average  for 
independent  symbols  drawn  from  the  same  constellation  as  the  training  sequence,  i.e., 
E{ak,iaJ^}  =  Ski-  Numerical  evaluation  shows  that  this  approximation  is  very 
accurate  even  for  training  sequences  as  small  as  L  =  10  or  20.  For  the  FH  model  in 
equation  35,  this  leads  to  the  following  approximations  that  are  used  in  the  CRB  and 
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ZZB  expressions  in  the  following  section.  First  we  define  the  autocorrelation  function 
of  the  (real-valued)  pulse  shape, 


hit  —  hit)  dt 


=  sine 


cos  (nat/T) 

1  -  {2at/Tf ' 


(49) 

(50) 


where  equation  50  is  the  raised-cosine  pulse  with  excess  bandwidth  (or  rolloff  factor) 
a  G  [0, 1],  which  is  the  result  when  hit)  is  a  square-root,  raised-cosine  pulse.  Note 
that  ph(0)  =  1.  Then  the  signal  autocorrelation  functions  in  equation  38  have  the 
form 


PiiO  =  LphiO  +  '^apittliPh  i^  +  il-  k)T) 

k^l 

-  LphiCj,  f  =  l,...,iV.  (51) 


Equation  47  has  the  form 
j  Siit)  Siit)*  dt  =  L  j 

f  =  l,...,iV, 


hit)  hit)  dt  +  i  /  hit  —  kT)  hit  —  IT)  dt 

k^i 


(52) 

(53) 


where  the  first  term  in  equation  52  is  zero  due  to  the  even  and  odd  symmetry  of  hit) 
and  hit),  respectively,  and  the  second  term  is  approximately  zero  due  to 
approximating  the  time  average  over  the  symbols  by  the  ensemble  average 
EWk,iCdii}  =  5ki-  The  signal  bandwidth  in  equation  48  has  the  form 


|s*(t)|^  dt  =  p(0)  (27ri?s 


=  L 


hit)  dt  +  cbk,icdi^i  /  hit  —  kT)  hit  —  IT)  dt 

ki^i 


L  i2T[Bhf  ,  i  = 


(54) 


where  Bh  Hz  is  the  mean-squared  bandwidth  of  the  pulse  hit).  Note  that  p(0)  =  L, 
so  the  signals  Sj(t)  have  the  same  bandwidth  as  the  pulse  hit),  as  expected  for  linear 
modulation.  When  hit)  is  the  square-root,  raised-cosine  pulse  with  excess  bandwidth 
(or  rolloff  factor)  a  G  [0, 1],  the  mean-squared  bandwidth  Bh  is  (i5) 


9  1  1 
TDZ  _ _ 

h-  t2  12 


1  + 


(55) 


Assumption  1  that  the  signal  energy  is  the  same  on  each  channel  in  equation  45  implies 
that  the  log-likelihood  (LL)  functions  for  the  conditional  and  unconditional  cases  are. 
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respectively, 


Lc  {R  I  7;  r) 
Cu  {R  I  r) 


1 

K 


5^[2y7“(r)-h..|y(0)] 


i=l 


_ fi _ y  y 

K  A/;  +  ajp(0)^  • 


NONCOH , 


(56) 

(57) 


These  LL  functions  are  used  in  the  next  section  to  derive  the  MLEs,  CRBs,  and  ZZBs  for 
the  various  channel  models. 


5.  TDE  Estimators  and  Bounds 

In  this  section,  the  MLEs,  CRBs,  and  ZZBs  are  presented  for  the  hve  channel  models 
summarized  in  table  1. 

5.1  Maximum  Likelihood  Estimators 

We  begin  with  the  MLEs,  and  we  dehne  the  sums  of  coherent  and  noncoherent  matched 
Liters  as  follows, 

=  5^  Re  (  /  O'-  s,{t  -  0*  nit)  dt]=Y.  (®) 

i=l  ^  i=l 

V  „  2  N 

^NONCOH(j)_^  s,{t-iYn{t)dt  (59) 

i=l  i=l 

For  case  1,  the  MLE  is  the  value  of  r  that  maximizes  the  LL  in  equation  56  when  7  is 
known.  The  second  term  in  equation  56  is  independent  of  r,  so  the  MLE  for  case  1  is  the 
value  that  maximizes  the  sum  of  coherent  matched  Liters, 

fi  =  arg  max  (60) 

For  case  2,  7  is  an  unknown,  deterministic  parameter,  so  the  MLE  requires  the  joint 
maximization  of  the  LL  in  equation  56  with  respect  to  r  and  7.  For  a  given  TDE  estimate 
r,  the  MLE  of  each  7*  is 

%  =  J  Siit-T)*ri{t)dt,  f  =  l,...,iV.  (61) 
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The  estimates  7*  in  equation  61  are  substituted  for  7  in  the  LL  of  equation  56,  yielding  the 
MLE  for  r  in  case  2  as  the  value  that  maximizes  the  sum  of  noncoherent  matched  filters, 

f2=  argmax 

For  case  3,  7  ~  CN  (0,  /)  so  the  MLE  for  r  is  the  value  that  maximizes  the 

unconditional  LL  in  equation  57,  which  is  identical  to  equation  62: 

f3=  argmax 

For  cases  lA  and  2A,  the  MLEs  are  given  by  fi  in  equation  60  and  f2  in  equation  62, 
respectively.  The  MLE  performance  is  different  in  cases  lA  and  2A  (compared  with  cases  1 
and  2  in  which  7  is  fixed)  due  to  averaging  over  realizations  of  values  for  7  ~  CN  (0,  I). 

The  MLEs  are  summarized  in  table  1,  where  “Coh  ME”  refers  to  equation  58  and 
“NonCoh  ME”  refers  to  equation  59. 

5.2  Cramer- Rao  Bounds 


Next  we  consider  the  CRB.  In  general,  let  us  denote  the  deterministic  parameters  by  the 
vector  ©  and  the  LL  function  by  £(R  |  ©).  Then  the  FIM,  J,  has  elements  for  parameters 
6*„  G  ©  given  by 

The  CRB  on  6n  is  the  nth  diagonal  element  in  the  inverse  of  the  FIM,  CRB(6*„)  =  [J~^]n,n- 


For  cases  1,  2,  and  3,  r*  is  a  complex,  Gaussian  random  vector,  r*  ~  CN(^j,  G*), 
independent  for  i  =  1, . . . ,  A^,  for  which  the  FIM  elements  are  given  by  {21) 


N 


i=l 


tr 


d9„ 


d  9  r. 


+  2  ■  Re 


djM 

d  9m  *  d9„ 


(65) 


where  tr(-)  denotes  the  trace  of  a  matrix.  For  cases  1  and  2,  /x*  =  7jSj(r)  and  C*  =  J\foI, 
with  ©  =  [r]  for  case  1  and  ©  =  [Re  (7) ,  Im  (7) ,  r]  for  case  2.  For  case  3,  fii  =  0  and 
Cj  =  Sj  in  (42),  with  ©  =  [r]. 


The  CRB  for  cases  1  and  2  can  be  obtained  from  the  following  FIM  for  case  2, 


/ 

Re  (7) 

\  2 

p{0)Inxn 

0  0 

J2 

Im  (7) 

j 

0 

p{0)Inxn  0 

V 

r 

0 

0  p(0)  (27rR.)^Ef=ilTp. 

(66) 
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where  assumptions  1,  2,  and  3  in  equations  45,  47,  and  48  have  been  used.  The 
block-diagonal  structure  of  equation  66  implies  that  the  parameters  7,  r  are  decoupled  in 
the  FIM,  so  the  CRBs  for  r  are  identical  for  cases  1  and  2, 


CRBi(r|7) 

CRB2(r|7) 


1  ^ _ 

['^2]2Ar+i^2V-|-l  2(2,rB,)7(0)  Eii  hip 

1  AC, 


p(0)Ef.i  l7iP' 


(67) 

(68) 


The  notation  on  the  left  side  of  equations  67  and  68  is  used  to  emphasize  that  the  CRBs 
depend  on  the  particular  values  of  the  channel  gains,  7.  The  FIM  for  case  3  is  obtained 
from  equation  65  with  =  0,  C*  =  S*  in  equation  42,  and  ©  =  [r].  Using  assumptions  1, 
2,  and  3  in  equations  45,  47,  and  48,  the  CRB  for  case  3  can  be  expressed  as 


CRB3  (r) 


1  +  p(0)  {all Mo) 

2  (27rR,)'  N  [p{Q)allNof' 


(69) 


The  CRB  in  equation  69  is  the  “true”  CRB  for  the  Rayleigh  fading  channel  model.  The 
derivation  of  equation  69  is  tractable,  but  it  requires  signihcantly  more  computation  than 
the  CRBs  in  equations  67  and  68  in  which  7  is  a  deterministic  parameter.  A  number  of 
bounds  that  are  looser  than  the  CRB  have  been  studied  for  models  with  random  nuisance 
parameters,  such  as  the  MCRB,  hybrid  CRB  (HCRB),  ACRB,  and  the  Extended 
Miller-Chang  Bound  (EMCB)  {13,  12,  37).  All  of  these  bounds  begin  with  a  FIM  in  which 
7  is  modeled  as  a  deterministic  parameter,  and  then  an  averaging  operation  is  performed 
over  7  on  some  function  of  the  FIM.  The  CRB  is  not  dehned  for  channel  models  lA  and 
2A  in  table  1  (case  3  provides  the  “true”  CRB  for  Rayleigh  fading),  but  the  MCRB  and 
ACRB  are  closely  related  to  cases  lA  and  2A,  respectively.  The  MCRB  is  similar  to  our 
channel  model  lA,  where  the  FIM  element  in  the  denominator  of  equation  67 
corresponding  to  known  7  is  averaged  {12),  yielding 


CRBi^  (r)  ^  MCRB  (f)  = 


'2j2Ar-|-l,2V-|-l 


2  (27rR,)'  p(0)iVa2 


The  ACRB  {37)  is  similar  to  our  channel  model  2A,  where  the  quantity  in  the  denominator 
of  equation  68  corresponding  to  unknown,  deterministic  7  is  averaged,  yielding 


CRB2U  (r)  ^  ACRB  (f)  = 


J2V] 


2V-|-1,2V-|-1 


Mo 

2  {2TiBsf  p(0)iVa2‘ 


It  is  shown  in  {12)  that  the  ACRB  is  a  lower  bound  on  the  EMCB  via  Jensen’s  inequality. 
Note  from  equations  69,  70,  and  71  that  the  MCRB  =  ACRB  are  looser  than  the  “true” 
CRB  for  Rayleigh  fading  in  equation  69,  and  that  the  three  bounds  coincide  at  high  SNR; 


CRBs  (r)  >  MCRB  (r)  =  ACRB  (f)  (72) 

.  f  — ^00 

CRBs  (r)  ^  MCRB  (r)  =  ACRB  (f) .  (73) 
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The  CRB  results  are  specialized  to  the  baseband  FH  signal  model  in  equation  35  as 
follows.  For  channel  model  cases  1  and  2,  the  (deterministic)  mean  SNR  per  symbol  over 
the  N  channels  is 

SNR*.  (7)  =  ^5:^.  (74) 

i=l  ° 

and  for  cases  lA,  2A,  and  3  (Rayleigh  fading),  the  mean  SNR  per  symbol  is  dehned  as 

SNR=Ft^  =  ^  =  B^{SNR*.(7)}.  (75) 

Then  using  p(0)  =  L  from  equations  45  and  51,  Bg  =  Bh  from  equation  54,  and  the  SNR 
per  symbol  dehnitions  in  equations  74  and  75,  the  CRBs  in  equations  67-71  are  as  follows 
for  the  FH  signal  model, 

CRBi  (f  I  7)  =  CRB2  (r  I  7)  =  [2  (27r  B^f  N-L-  SNR^et  (7)]  (76) 

CRB3(r)=  2  {271  Bhf  N  (77) 

CRBi^  (r)  MCRB  (f)  =  [2  (27r  Bhf  N-L-  WR\  (78) 

CRB2A  (r)  ^  ACRB  (r)  =  [2  (27r  B^f  N-L-  WR\  .  (79) 

We  observe  the  following  from  the  CRB  analysis. 

•  The  time  delay  r  and  channel  gains  7  are  decoupled  in  the  FIM,  so  the  CRBs  are 
identical  in  cases  1  and  2  where  the  channel  gains  are  modeled  as  known  and 
unknown  deterministic  parameters,  respectively. 

•  The  MCRB  and  ACRB  for  channel  model  cases  lA  and  2A  are  identical  to 
CRBi  =  CRB2  without  fading,  as  long  as  SNRdet  =  SNR.  Furthermore,  the  true 
CRB3  for  Rayleigh  fading  coincides  with  the  other  CRBs  at  high  SNR.  Therefore  the 
CRB  is  not  very  sensitive  to  the  model  for  the  channel  gain  parameters,  particularly 
at  high  SNR. 

•  The  CRBs  are  inversely  proportional  to  the  SNR,  the  pulse  bandwidth  Bh,  and  the 
total  number  of  symbols,  Q  =  L  -  N,  where  N  is  the  number  of  hops  and  L  is  the 
number  of  symbols  per  hop.  The  only  exception  is  CRB3  in  equation  77,  which  for 
fixed  Q  =  LN  and  low  SNR  is  minimized  by  A^  =  1  hop.  The  other  CRBs,  as  well  as 
CRB3  at  high  SNR,  depend  only  on  the  total  number  of  symbols,  so  the  CRB  does 
not  reveal  a  diversity  advantage  when  more  hops  are  used  with  fewer  symbols  per  hop. 

We  will  show  in  the  following  subsection  that  the  Ziv-Zakai  bounds  are  quite  sensitive  to 
the  channel  model,  and  that  when  the  channel  gains  are  random  variables  (fading),  then 
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greater  diversity  (larger  N)  provides  significant  improvement  in  the  threshold  SNR  at 
which  the  MLE  performance  is  equal  to  the  CRB.  The  ZZBs  indicate  that  with  fading 
channels,  the  CRB  on  TDE  is  a  tight  bound  only  when  there  is  sufficient  diversity  {N  >  1) 
and  the  SNR  exceeds  a  threshold  value. 


5.3  Ziv-Zakai  Bounds 


The  four  distinct  ZZBs  summarized  in  table  1  are  developed  in  this  section  using  the 
expressions  in  equations  13-17.  The  general  hypothesis  test  in  equation  13  is  stated  as 
follows  for  the  model  of  TDE  over  N  parallel  channels  in  equation  34: 

Ho  :  ri{t)  = -  a)  +  i  =  1, . . . ,  N  (80) 

Hi:  ri{t)  = 'yiSi{t  -  {a  +  9))  +  ni{t),  i  =  (81) 

The  hypotheses  are  modeled  as  equally  likely.  In  the  ZZB  formulation,  r  is  modeled  as  a 
uniform  random  variable  on  the  interval  [—D/2,  D/2],  which  represents  a  priori  knowledge 
about  the  range  of  possible  values  for  the  time  delay  parameter,  r.  Therefore  in  equations 
80  and  81,  6*  >  0  and  a,  {a  +  9)  G  [—D/2,  D/2].  Let  Pe{a,  a  +  9)  denote  the  minimum 
probability  of  error  in  deciding  between  Hq  and  Hi  for  particular  values  of  a,  a  +  9.  The 
Pe{a,  a  +  9)  is  derived  from  a  likelihood  ratio  test  (LRT),  so  it  depends  on  the  model  for 
the  channel  gain  vector,  7.  The  conditional  LL  in  equation  56  is  used  to  find  Pe  for  channel 
models  1  and  lA,  with  7  fixed  in  case  1  and  averaged  over  7  in  case  lA.  For  case  2,  a 
generalized  LRT  (GLRT)  formed  with  ML  estimates  of  the  channel  gain  parameters  in 
equation  61  substituted  into  equation  56  reduces  to  a  LRT  based  on  equation  57,  with  7 
fixed.  Case  2A  includes  subsequent  averaging  over  7.  The  Pg  for  case  3  is  based  on  the 
unconditional  LL  in  equation  57,  and  is  found  to  be  identical  to  case  2A.  This  summarizes 
the  procedure  used  to  derive  the  four  ZZBs  in  table  1,  and  the  details  are  presented  next. 


For  case  1,  an  error  occurs  in  the  LRT  based  on  equation  56  in  the  following  cases: 


N 

N 

<^yCOH, 

[a  -\-  9) 

(82) 

i=l 

i=l 

N 

N 

i=l 

>^yCOH, 

i=l 

[a  -\-  9) 

(83) 

where  is  defined  in  equation  39.  The  minimum  probability  of  error  for  case  1, 

Pe,i{(i,  0,  +  9  \  ^),  is  one-half  the  probability  of  equation  82  plus  one-half  the  probability  of 
equation  83.  It  is  straightforward  to  show  that  equations  82  and  83  have  equal  probability 
that  is  independent  of  a,  so  for  case  1  with  7  known  and  a  =  0, 

^yCOH(o)  _yCOH(^)j 


N 


Pg,i  (0|7)=Pr<^X] 


i=l 
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The  result  for  case  lA  is  obtained  by  averaging  over  7, 


(«)  =  £:{P.J  («|7)}. 


(85) 


For  case  2,  the  hypothesis  test  equations  80  and  81  contains  the  unknown,  deterministic 
channel  gain  parameters,  7.  The  GLRT  is  formed  by  substituting  the  ML  estimates  of  the 
channel  gain  parameters  in  equation  61  into  equation  56,  which  by  equation  62  results  in 
LL  functions  that  are  a  sum  of  noncoherent  matched  Liters,  as  in  equation  57.  As  in  case  1, 
the  error  has  equal  probability  under  Hq  and  Hi  and  is  independent  of  a,  so  with  a  =  0 

^yNONCOH^Q^2  _yNONCOH^^^2j  |  |  ^ 

where  jg  fjegj^ed  in  equation  40.  The  result  for  case  2A  is  obtained  by 

averaging  over  7, 


N 


(0|7)=Pr  X] 


i=l 


Pe,2A  {9)=E{Pe,2  (^  |  7)  }  • 


(87) 


For  case  3,  the  LRT  is  based  on  equation  57,  the  error  has  equal  probability  under  Hq  and 
Hi  and  is  independent  of  a,  so  with  a  =  0 

Pe,3  {9)  =  Pr  |x]  -  yNONCOH^^^2j  ^  g  |  |  ^gg^ 

=  Pe,2A  (9)  .  (89) 

The  probability  in  equation  88  is  computed  with  respect  to  the  additive  noise  and 
7  ~  CN  (0,  (Tg  I),  so  the  result  is  identical  to  case  2A,  as  indicated  by  equation  89.  Next, 
we  provide  expressions  for  the  probabilities  in  equations  84-86,  and  88.  The  ZZB  for  each 
case  is  obtained  by  using  the  corresponding  Pg  in  equation  17,  which  is  repeated  here  for 
convenience, 

ZZB  (r)  >  ^  /  9  V[{D  -  9)P49)]  d9,  (90) 

P  Jo 

where  V[-]  is  a  nonincreasing  function  that  fills  the  valleys  of  the  bracketed  function. 
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We  begin  with  case  1  and  the  evaluation  of  equation  84,  where 

U,  =  I  hfo  =  Re  [  7;  s,{tr  [7,  s,{t)  +  n,{t)]  dt 


=  l7*Pp(0)  +  Re 


li  /  Si(t)*ni(t)dt 


'X' 


V,  =  \Ho  =  Re  /  7*  s,{t  -  9)*  [7,  s,{t)  +  n,{t)]  dt 


=  |7*rRe[p*(0)]  +  Re 


li  /  Si{t  -  9)*  ni{t)  dt 


~b' 


(91) 

(92) 

(93) 

(94) 


The  Ai,  Bi  are  complex,  circular,  correlated,  Gaussian  random  variables  with  zero  mean 
and 


E{\AX}=p{0)Afo 

e{\b/}=p{o)K 
E{AB*}  =  p,{-9)K 

1  p*(^)7p(o) 

p*7)/p(0)  1 


B, 


CN  0,  piO)Afo 


It  follows  that  the  second  terms  in  equations  92  and  94  are  real-valued,  correlated, 
Gaussian  random  variables  with  zero  mean  and  covariance  matrix 


Re  bi  dli] 
Re  b*i  Bi] 


~  N 


1 

2 


\Y?P{^)^o 


1 

Re[p.(P)]/p(0) 


Re[p.(P)]/p(0) 

1 


(95) 

(96) 

(97) 

(98) 


(99) 


and  the  random  vectors  in  equation  99  are  independent  for  i  =  1,  2, . . . ,  iV.  The  minimum 
probability  of  error  for  case  1  in  equation  84  is  then 


Pe,l  b  I  7) 


N 


Pr\J2[U.-Vi\ 


<  0 


i=l 


(100) 


Using  equations  92,  94,  and  99,  the  independence  over  i,  and  the  assumption  Pib)  =  p(7 
in  equation  46, 

N  /  N  N  \ 

Y^[U.-Vi\  ~N  (p(0)-Re[p(9)])  (p(0)  -  Re  [p(9)])  AC,  ^  ,  (101) 

i=l  V  *=1  *=1  / 

SO  equation  100  is  expressed  as 

Pe,i  {0\A=Q  (7p(0)  ■  N  ■  SNRdet  (7)  ■  (1  -  Re  W)]  /p(0)))  (102) 
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where 


Q{y)  = 


exp(— m^/2)  du 


(103) 


and  SNRdet  (7)  is  defined  in  equation  74.  The  probability  of  error  for  case  lA  in  equation 
85  can  be  evaluated  using  {60), 


n{e) 


PeM  («) 


.^^p(0)-SNR-(l-ReH«)|/p(0)). 
(l-pW)"g(/V-l  +  i) 


(104) 

(105) 


where  SNR  is  defined  in  equation  75. 


Next,  we  follow  similar  reasoning  for  the  LRT  based  on  the  sum  of  noncoherent  matched 
filters  to  evaluate  equations  86  and  88.  We  define 


Wi  =  j  Si{t)*  ['yiSi{t)  +  ni{t)]  dt  (106) 

=  7p(0)+A  (107) 

y)NONCOH^Q^2  I  _  |^,|2 

ATi  =  J  Si{t  -  9)*  ['jiSi{t)  +  ni{t)]  dt  (109) 

=  liPi{9)  +  Bi  =  'jip{9)  +  Bi  (110) 

yNONCOH('^^2  I  ^  |^,|2  ^ 


where  Ai,  Bi  are  defined  in  equations  92  and  94  and  are  distributed  as  in  equation  98.  We 
have  assumed  pi{^)  =  p{^)  as  in  equation  46  for  the  rightmost  equality  in  equation  110. 

The  probability  of  error  for  case  2  in  equation  86  is  then  evaluated  as 

Pe,2  (^|7)  =Pr 

where  7  is  a  deterministic  parameter.  For  cases  2A  and  3,  the  probabilities  in  equations  87 
and  88  are  evaluated  as 

W|'-|W|2]  <o|  (113) 

where  7  ~  CN  (0,(T^/).  In  equations  112  and  113,  the  random  variable 

0^*1^  ~  is  ^  sum  of  quadratic  forms  in  which  the  N  pairs  {Wi,Xi)  are 

complex,  correlated,  Gaussian  random  variables  that  are  independent  for  i  =  1, . . . ,  N.  For 


Pe,2A  {0)  =  Be, 3  {9)  =  Pr 


N 


[iw-y  -  ivi"] 


<  0 


(112) 


i=l 
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case  2  in  which  7  is  deterministic,  the  distribution  of  {Wi,Xi)  follows  from  equations  107, 
110  and  98  as 


W, 

X, 


~  CN 


1 

p(^)/p(0) 


,  p(0)K 


1  p(0r/p(o) 

p(P)/p(0)  1 


For  case  3  in  which  7  CN  (0,(T^/),  the  distribution  of  (lFj,Xj)  is 


W, 

X, 


CN  0,  p(0)  K  [1  +  SNR  ■  p(0)' 


1 

p(^)/p(o) 


m*/p(o) 

1+SNR.(|p(0)|Vp(O)) 

1+MR-p(0) 


(114) 


(115) 


The  probabilities  in  equations  112  and  113  can  both  be  evaluated  using  the  general  result 
in  (59)  or  (60),  which  is  applicable  with  nonzero  mean  as  in  equation  114  as  well  as  zero 
mean  in  equation  115.  Both  cases  can  be  unihed  by  a  Rician  fading  model,  in  which  case 
Pe  can  be  obtained  from  (59),  (60),  and  then  equations  112  and  113  would  be  special  cases. 
It  would  be  useful  to  work  out  the  Rician  result,  but  here  we  restrict  attention  to  Rayleigh 
fading  (case  3),  and  equation  113  has  the  form  (59),  (60) 


PeS  {«) 


4(l  +  SNR-p(0))  " 

_  (p(O).SNR)".  [l-|p(»)/p(0)|2], 

/I -!/(«) y 2iV-l  \  p  +  i^(9)y 

I  2  j 


(116) 


(117) 


We  conclude  this  section  by  specializing  the  ZZB  results  to  the  FH  signal  model  in 
equation  35,  where  p(0)  =  L  from  equation  45  and  p(^)  ~  L ph(i)  from  equation  51.  The 
ZZB  for  each  case  is  numerically  evaluated  from  equation  90  using  the  following  Pg 
expressions,  where  the  autocorrelation  of  the  pulse  shape,  Ph(i))  is  real- valued. 


Pep  (e  1 7) 

p(e) 


Pe,lA  (0) 


Pe,2A  (9)  =  Pe,3  (0) 


q[^L-N-  SNRdet  (7)  ■  (1  -  Ph(e))) 

1  1 

L -MR- (I- ph(9))_ 

^  4  (1  +  L  ■  Mr)  ’ 

(L-MR)^[1-p,(P)2]_ 

(i-v(e)V^-^'^s^  (  2N  -1\  (i  +  v(e)V 

V  2  J 


(118) 

(119) 

(120) 
(121) 
(122) 


These  Pg  expressions  can  be  used  to  analytically  approximate  the  SNR  thresholds  for  TDE 
on  parallel,  fading  channels.  In  the  next  section,  simulation  results  are  presented  that 
compare  the  MSE  performance  of  the  MLE  with  the  CRBs  and  ZZBs  for  various  cases. 
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6.  Examples  and  Simulation  Results 


This  section  presents  numerical  evaluation  of  the  bounds  (CRBs  and  ZZBs)  for  various 
cases,  and  we  also  compare  the  MSE  performance  of  the  MLE  with  the  bounds  for  various 
scenarios.  We  consider  the  baseband  EH  signals  as  in  equation  35  with  symbol  time  is 
Ts  =  10“^  s  and  pulse  shape,  h(t),  a  square-root,  raised-cosine  with  0  excess  bandwidth 
(a  =  0).  The  known  symbols  are  generated  as  iid  QPSK  sequences,  and  the  symbol 
sequence  on  each  hop  is  kept  fixed  over  all  simulation  runs  to  match  the  model  used  in  the 
analysis. 

We  begin  by  considering  a  single  symbol  per  hop,  so  =  0  and  L  =  1  in  equation  35,  which 
corresponds  to  TDE  in  equation  34  with  signals  Si{t)  =  h{t)  over  N  parallel  channels.  We 
evaluate  the  bounds  and  MLE  performance  as  N  is  increased,  and  we  show  that  the  ZZB  is 
a  tight  bound  for  the  MSE  of  the  MLE.  Note  that  for  this  case,  the  signal  autocorrelation 
in  equation  51  is  exact.  Next,  we  consider  the  EH  model  in  equation  35  with  the  total 
number  of  symbols  hxed,  so  Q  =  LN  is  constant.  The  CRB  depends  primarily  on  Q,  but 
we  show  that  the  ZZB  is  strongly  dependent  on  the  number  of  diversity  channels,  N,  and 
that  the  CRB  is  not  tight  when  N  is  small  or  when  the  SNR  is  less  than  a  threshold  value. 

6.1  One  Symbol  Per  Hop 

Figure  2  shows  the  simulated  MSE  in  comparison  with  the  CRB  and  ZZB  for  iL  =  0  and 
L  =  1  symbol  per  hop.  In  the  legend  of  figure  2(c),  “Coh  &  No  Fade”  corresponds  to 
channel  model  case  1,  “Coh  &  Rayleigh”  corresponds  to  case  lA,  “NonCoh  &  Rayleigh” 
corresponds  to  cases  2A  and  3  for  the  ZZB,  and  “Rayleigh”  corresponds  to  case  3  for  the 
CRB.  The  MSE  performance  is  shown  for  the  MLEs  in  cases  1,  lA,  and  3.  Note  in  part  (a) 
that  with  fading  (the  black  and  red  lines)  and  no  diversity  (A^  =  1),  the  ZZB  is 
significantly  larger  than  the  CRB,  and  the  CRB  is  not  achievable  an  any  SNR.  Parts  (b) 
and  (c)  show  that  with  diversity  {N  =  3  and  5),  the  ZZB  converges  to  the  CRB  when  the 
SNR  is  exceeds  a  threshold,  and  the  threshold  of  the  ZZB  is  close  to  the  threshold  of  the 
MSE,  particularly  for  the  cases  of  no  fading  (blue  line)  and  Rayleigh  fading  (black  line). 
The  MSEs  are  computed  based  on  2,000  Monte  Carlo  runs. 


TDE  ERROR  (sec) 


SQUARE  ROOT  OF  ZZB,  CRB,  &  SIMULATED  MSE:  N  =  1  CHANNELS 


SQUARE  ROOT  OF  ZZB,  CRB,  &  SIMULATED  MSE:  N  =  3  CHANNELS 


(b) 


SQUARE  ROOT  OF  ZZB,  CRB,  &  SIMULATED  MSE:  N  =  5  CHANNELS 


Figure  2.  Comparison  of  CRB,  ZZB,  and  MSE  for  case  of  L  =  1  symbol  per  hop  for  =  1,  3,  and  5 
channels. 

Figure  3  shows  CRBs  and  ZZBs  for  the  same  model  with  N  =  1,2, ...  ,8,  except  that  part 
(a)  corresponds  to  coherent  matched  filters  with  no  fading  (case  1),  (b)  corresponds  to  the 
average  of  part  (a)  over  Rayleigh  channels  (case  lA),  and  (c)  corresponds  to  noncoherent 
matched  filters  and  Rayleigh  fading  (case  3).  Note  that  the  SNR  thresholds  differ 
considerably  in  parts  (a),  (b),  and  (c).  The  ZZB  converges  to  the  CRB  in  parts  (a)  and  (c). 
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TDE  ERROR  (sec) 


In  part  (b),  the  CRB  is  not  defined,  but  the  MCRB  and  ACRB  are  plotted.  Since  there  is 
no  “true”  CRB  for  part  (b),  it  is  not  surprising  that  the  ZZB  does  not  converge  to  the 
MCRB  =  ACRB. 


SQUARE  ROOT  OF  ZZB  &  CRB:  Coherent  &  No  Fading  SQUARE  ROOT  OF  ZZB  &  CRB;  Coherent  &  Rayleigh  Fading 


(a)  (b) 


SQUARE  ROOT  OF  ZZB  &  CRB:  Noncoherent  &  Rayleigh  Fading 


Figure  3.  CRB  and  ZZB  for  case  of  L  =  1  symbol  per  hop  for  =  1,  2, . . . ,  8.  (a)  Known  channel 
with  fixed  gain;  (b)  average  of  case  (a)  over  Rayleigh  channels,  so  the  channel  gain  is  known 
on  each  realization  and  used  in  a  coherent  matched  filter;  and  (c)  Rayleigh  channels  with 
noncoherent  matched  filters. 
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6.2  Multiple  Symbols  Per  Hop 


Figure  4  shows  the  CRB  and  ZZB  for  hxed  total  symbols  Q  =  LN  =  80,  for  N  =  1,  2, 4,  8 
channels  with  L  =  80, 40,  20, 10  symbols  per  hop,  respectively.  Note  that  the  CRB  depends 
primarily  on  Q  and  varies  with  N  only  at  low  SNRs  for  the  Rayleigh  fading  case  in  part 
(c).  In  part  (a),  with  no  fading  and  coherent  matched  hlters,  the  ZZB  and  CRB  depend 
only  on  Q  and  are  identical  as  N  is  varied.  In  part  (c),  the  ZZB  changes  signihcantly  as  N 
is  varied.  For  iV  =  1  with  no  diversity,  the  ZZB  is  much  larger  than  the  CRB.  For 
TV  =  2, 4,  8  in  part  (c),  the  threshold  SNR  decreases  as  the  number  of  hops  is  greater. 
However,  the  difference  between  iV  =  4  and  iV  =  8  is  small,  so  the  greatest  diversity  gain 
occurs  with  the  hrst  few  additional  channels. 

Figures  5  and  6  contains  simulation  results  based  on  4,000  runs  and  6,000  runs, 
respectively,  to  estimate  the  MSE  using  the  MLE  for  channel  model  lA  in  part  (a)  and 
channel  model  3  in  part  (b).  Part  (b)  is  the  case  of  practical  interest  with  noncoherent 
matched  hlters  and  Rayleigh  fading,  and  note  that  the  MSE  closely  tracks  the  ZZB  for 
=  4,  and  the  MSE  is  larger  than  the  ZZB  for  =  1  and  N  =  2.  In  all  cases,  the  ZZB  is 
a  much  tighter  bound  than  the  CRB. 
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TDE  ERROR  (sec) 


(a)  (b) 


SNR(dB) 

(c) 

Figure  4.  CRB  and  ZZB  for  fixed  total  number  of  symbols  Q  =  LN  =  80,  for  N  =  1,  2,  4,  8  so 
L  =  80, 40,  20, 10  to  show  the  tradeoff  of  shorter  dwells  per  hop  and  more  hops,  (a)  Known 
channel  with  fixed  gain  (channel  model  1);  (b)  average  of  case  (a)  over  Rayleigh  channels,  so 
the  channel  gain  is  known  on  each  realization  and  used  in  a  coherent  matched  filter  (channel 
model  lA);  and  (c)  Rayleigh  channels  with  noncoherent  matched  filters  (channel  model  3). 
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TDE  ERROR  (sec)  TDE  ERROR  (sec) 


(a) 


(b) 


Figure  5.  Comparison  of  CRB,  ZZB,  and  MSE  for  Q  =  LN  approximately  fixed,  with  =  1,  2, 4  and 
L  =  81,41,21.  (a)  Coherent  matched  filters,  averaged  over  Rayleigh  channel  realizations 
(channel  model  lA).  (b)  Noncoherent  matched  filters  with  Rayleigh  fading  (channel  model 
3).  MSEs  are  estimated  from  4,000  runs. 


(b) 


Figure  6.  Comparison  of  CRB,  ZZB,  and  MSE  for  Q  =  LN  approximately  fixed,  with  =  1,  2, 4  and 
L  =  81,41,21.  (a)  Coherent  matched  filters,  averaged  over  Rayleigh  channel  realizations, 
(b)  Noncoherent  matched  filters  with  Rayleigh  fading.  MSEs  are  estimated  from  6,000  runs. 


33 


7.  Concluding  Remarks 


Several  items  for  continued  research  are  suggested  by  this  work.  First,  it  would  be  useful  to 
generalize  the  fading  model  to  the  Rician  case,  and  then  the  fixed  and  Rayleigh  7  models 
would  be  special  cases  of  the  Rician  model.  The  extension  to  Rician  fading  is 
straightforward  based  on  the  CRB  and  ZZB  analysis  developed  in  this  report.  Second,  it 
would  be  useful  to  analytically  characterize  the  threshold  behavior  of  the  ZZB  with  respect 
to  SNR  and  the  number  of  diversity  channels,  N.  A  simple  formula  would  be  desirable  to 
complement  the  simulation  results.  Third,  the  baseband  FH  signal  model  was  investigated 
in  detail,  and  it  will  be  interesting  to  develop  corresponding  results  for  the  sum-of-tones 
model  in  equation  37. 
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ACRE 

asymptotic  CRB 

AWGN 

additive  white  Gaussian 

BE 

Barankin  bound 

CRB 

Cramer-Rao  bounds 

CSD 

cross-spectral  density 

EMCB 

Extended  Miller-Chang  Bound 

FH 

frequency-hopping 

FIM 

Fisher  Information  Matrix 

GCC 

generalized  cross-correlator 

GIRT 

Generalized  LRT 

HCRB 

hybrid  CRB 

LL 

log-likelihood 

LRT 

likelihood  ratio  test 

MCRB 

modified  CRB 

MFs 

matched  filters 

MLE 

maximum-likelihood  estimator 

MSE 

mean-squared  error 

pdf 

probability  density  function 

PSD 

power  spectral  density 

RMS 

root-mean-square 

SNR 

signal-to-noise  ratio 
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TB 

time-bandwidth 

TDE 

time-delay  estimation 

TDOA 

time  difference  of  arrival 

WWBs 

Weiss-Weinstein  bounds 

ZZB 

Ziv-Zakai  bound 
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